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ABSTRACT 

Roughly perpendicular to SS433's famous precessing jets is an outflowing "ruff" of radio-emitting 
plasma, revealed by direct imaging on milli-arcsecond scales. Over the last decade, images of the 
ruff reveal that its orientation changes over time with respect to a fi xed sky co-ordinate grid. F or 



example, during two months of daily observations with the VLB A by Mioduszewski et al. (2004), a 
steady rotation through ~ 10° is observed whilst the jet angle changes by ~ 20°. 
The ruff reorientation is not coupled with the well-known precession of SS433's radio jets, as the ruff 
orientation varies across a range of 69° whilst the jet angle varies across 40°, and on greatly differing 
and non-commensurate timescales. 

It has been proposed t hat the ruff is fed by SS433s circumbinary disk, discovered by a sequence of 



optical spectroscopy by Blundell et al. (2008), and so we present the results of 3D numerical simulations 
of circumbinary orbits. These simulations show precession in the longitude of the ascending node of all 
inclined circumbinary orbits — an effect which would be manifested as the observed ruff reorientation. 
Matching the rate of ruff precession is possible if circumbinary components are sufficiently close to 
the binary system, but only if the binary mass fraction is close to equality and the binary eccentricity 
is non-zero. 

Subject headings: Stars: Binaries: Close, Stars: Individual: SS433 



1. INTRODUCTION 

Imaging on milli-arcsecond scales is routine only at ra- 
dio wavelengths with very long baseline interferometry 
(VLBI) techniques. Such techniques, applied to observa- 
tions of the microquasar SS433 reveal extended emission 
roughly perpendicula r to its famous precessing, knotty 
100 AU dBlundell et al.|[2001| |Mioduszewski 



et al. 


20041 


SS43a 


s ruff. 



A succession of optical spectra of the Balmer H-a line 
revealed a p air of stationary com ponents that were in- 



terpreted by Blundell et al. (2008) as arising from a cir 
cumbinary disk (because of their lack of depen dence on 



the orbital m otion of the binary components). Blundell 



et al.| (120081) further posit that this circumbinary disk 
feeds the radio ruff. 

We have performed fully 3-D numerical simulations of 
circumbinary test particles in orbit around synthetic bi- 
nary systems as a function of mass fraction and orbital 
eccentricity. Our simulations reveal clear trends in the 
behaviour of the plane of a circumbinary orbit with re- 
spect to the plane of the binary which we compare with 
time-resolved observations of SS433's ruff. 

2. THE CHANGING ORIENTATION OF SS433'S RUFF 

We extracted an estimate of the mean angle of the over - 
all ruff from images presented by [Blundell et al. (2 001|); 



al. 



^ fl 2QQ2a N l 1999 l > and |Miodusze^i"eT^ 
covering VLB! observations of SS433 from 1995 
to 2003. We show in Figure [I] a subset of these images 
that clearly demonstrates the change in angle, projected 
on the plane of the sky, of SS433's ruff with time. 

From each of the VLBI images available to us we ex- 
tracted this "ruff angle", measured anti-clockwise from 
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Fig. 1. — : Examples of radio images on milli-arcsec 
scales showing the famous precessing jets of SS433 and the 
differently precessin g ruff - highlighted in green. 
Top: 7 March 1998, |Blundell et al.| f2001) (their fig. 2); 
Centre: 6 June 1998, Paragi et al. (2002a) (their fig 3); 
Lower: 13 Feb 2000, |Paragi et al.| (|2u02b) (their fig 3a). 
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celestial north. Where more than one published image 
was available for a given observation we measured all 
images. Higher resolution images reveal inner (more re- 
cently launched) ruff material, compared to the more ex- 
tended ruff emission visible in lower-resolution images. 

We infer the epoch of emission (hereafter referred to as 
the "corrected" date) via the extent of the features of the 
observable ruff assuming a distance to SS433 of 5.5 kpc 
(Blu ndell fc Bowler|2Q Q4) and the average wind speed of 
10, 000 km/ s as reported by |Mioduszewski et al.| ([2004 ) , 
although the launch mechanism lor these observed high 
speeds remains to be established. These corrections are 
of order 30 days. 

In Fig [2] we plot the ruff angle as a function of corrected 
Julian Date: these data clearly show systematic varia- 
tions through a range of at least 50°. A fairly smooth, 
steady change is even revealed within the Mioduszewski 
movie observations of summer 2003. 

We also show in Figure [2] the best fitting sinusoid. This 
fit is only to the Mioduszewski data-points, giving an ini- 
tial guess for the amplitude of 35° with no vertical offset, 
after which frequency, phase, amplitude and vertical off- 
set are fitted. 

It is remarkable that many of the earlier data points 
lie (without any further fitting) on this sinusoid, which 
shows a periodicity in the apparent angle of the ruff of 
552.5 days, or 42.2 binary orbital periods 2 . 



Paragi May 1995 
Blundell March 1998 
Paragi June 1998 
Paragi Feb 2000 
Mioduszewski Movie Obs 2003 




Fortunately, the time sampling of the available data 
enable us to rule out such correlations. In Fig|3]we re- 
plot the measured ruff angles folded over SS433's three 
major periodicities and show an absence of a preferred 
ruff orientation with orbital, nodding or precession phase. 
We note that the periodicity plotted in Fig[2] of 552.5 
days is neither commensurate with the precession period 
of 162.375 days, orbital period of 13.08 days, nor nodding 
period of 6.06 days. 

We also remark that the fit to the ruff angle varies 
through 69°, whereas the jet angle projected on the plane 
of the sky only has a range of 40° (Hjellming & Johnston 
1981). 



Precession Phase 



Orbital Phase 



Nodding Phase 



Fig. 3. — : Plot of observed ruff angle versus precession, 
orbital and nodding phases: there is no evidence of any pre- 
ferred ruff angle with any of these phases. Legend as Fig. [2] 



Fig. 2. — : The data points are ruff angle vs corrected Ju- 
lian Date (see §|2j), measured from all VLBI images of SS433 
a vailable to us. The sinu soidal line was fitted only to the 



4. DESCRIPTION OF THE 3-D SIMULATIONS 



Mioduszewski et al. (2004) data; details in text 



3. 



INDEPENDENCE OF RUFF MOTION FROM OTHER 
PERIODICITIES 



We first consider whether the apparent precession of 
SS433's ruff simply reflects other known periodicities in 
the SS433 system. There are three principal periodicities 
to consider: SS433's jet axis tra ces out a cone of semi- 
angle - 20° every 162.375 days flEikenberry et al ||2001[ ) 
known as its precess ion period; the orb ital period of the 
binary is 13.08 days ( |Kemp et al.(l 986); and third, there 
is a nodding superimposed on the precession of the jet 
axis with a period of 6.06 days, manifested a s beats with 
the p recession period at 5.83 and 6.28 days ( |Katz et al.| 
1982). 



One orbital period = 13.08 days ( Kemp et al.|l986| 



Blundell etH] fl2008D posit that SS433's ruff is fed 



from its circumbinary disk, and we further posit that 
the change in ruff angle arises because of an evolution in 
the orientation of this circumbinary structure. This mo- 
tivated us to investigate the behaviour of circumbinary 

orbits. 

We follow on from the 2D simulations of Holm an~fcl 
Wiegert ( 1999 ) by numerically integrating a suite of non- 



interacting test particles around a binary system of arbi- 
trary eccentricity and mass ratio, in fully 3D calculations. 

We solve the orbit of the binary directly by applying 
a dual bisection and Newton-Raphson algorithm to Ke- 
pler's equation. We apply an adaptive step-size fourth 
and fifth-order Runge-Kutta integrator to integrate the 
system for 10 4 binary orbital periods. During integration 
we remove test particles which come too close to either 
star, or escape the system. 

A general orbit in 3D about a centre of mass 
may be described by the Keplerian orbital elements 



SS433's orbital eccentricity 
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(a, e, z, W, w, v) as illustrated by Figured The eccen- 
tricity e, semi-major axis a and true anomaly v describe 
the motion of a test particle in its orbital plane, whilst the 
inclination z, longitude of the ascending node W and ar- 
gument of perihelion w describe the orientation of the or- 
bital plane of the test particles with respect to the plane 
of the binary. 



Figure [5] is a plot of the longitude of the ascending node 
and "measured angle" of an example test particle's orbit, 
as a function of binary eccentricity, binary mass fraction 4 
(alpha) and time in units of binary period (T). This 
particle orbits at radius 3.0 times the semi-major axis of 
the binary system and has inclination angle i = 36° to 
the plane of the binary. 




Fig. 4. — : A representation of the Keplerian orbital elements 

We track the orbital elements of each test particle 
about the centre of mass of the binary, outputting snap- 
shots as a function of time to a MySQL database 3 . We 
deem an orbit to be stable if it remains circular, but we 
notice that all inclined orbits precess in the longitude o f 
the ascending node in a kozai-like manner ( |Kozai|[l962| ); 
like a coin spinning on a tabletop (where the face of the 
coin is the plane of the circumbinary orbit). 

5. FINDINGS FROM OUR SIMULATIONS 

Since SS433 is an eclipsing system we know that the 
orbit al pla ne lies nearly perpendicular to the plane of the 
sky (£6.1), so we set the plane of our simulated binary 
thus. We project the orbit of a test particle of inclination 
i and longitude of the ascending node W onto the plane 
of the sky and extract the angle that the orbit appears to 
subtend relative to the edge-on orbital plane. Hereafter 
this is referred to as the "measured angle" . 

If SS433's ruff is d riven from circumbinary material 
( Blundell et al.|[2QQ8[ ) then the changing angle presented 
by the ruff is manifested in our simulations as this "mea- 
sured angle". As a test particle's orbital elements evolve 
over time it is the "measured angle" of its orbit that we 
contend is the ruff angle. 

3 These snapshots are of "osculating" orbits: the Keplerian or- 
bit about a central body that a test particle would have if other 
perturbations were not present. 




Fig. 5. — : Tracking test particles: we plot the longitude of 
the ascending node and "measured angle" of a circumbinary 
orbit as a function of time, for example values of (a) binary 
mass fraction, and (b) binary eccentricity. 

A full description of the parameter space explored by 
our simulations and the diagnostic machinery this pro- 
vides for the interpretation of circumbinary orbits will be 
presented in a forthcoming paper. Here, we focus on the 
trends pertinent to the changing ruff angle in SS433: 

• The closer the binary mass fraction to equality, the 
shorter the period of oscillation of the measured an- 
gle, i.e. precession of the longitude of the ascending 
node W. 

• The greater the binary eccentricity, the shorter the 
period of oscillation of the measured angle. 

• Fig|5]3 also shows that increasing binary eccentric- 
ity produces a deviation away from a steady pre- 
cession of the longitude of the ascending node W, 
and hence a departure from a pure sinusoidal curve 
in the oscillation of the measured angle. It is possi- 
ble that the very small fraction of data points that 
do not lie on the sinusoidal fit of Fig. [2] are deviant 
because SS433's binary orbit is eccentric. 

• The closer the test particle orbits to the binary 
system, the shorter the period of oscillation of the 
measured angle. This of course arises because of 
the stronger coupling of angular-momentum ex- 
change between closer-in particles and the binary, 
compared to particles that are further away. 

4 Binary mass fraction = Mi/(Mi+M2), where Mi is the lighter 
component and M2 the heavier binary component. 
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5.1. Determination of orientations within SS433 

From the fit shown in figure [2] we deduce that the ruff 
precesses with an amplitude of 34.5°. At the maxima 
and minima of the curve the plane of the circumbinary 
disk contains our line of sight, actually perpendicular to 
the plane of the sky. Since the orbital plane is also into 
the plane of the sky we determine that the circumbinary 
disk is therefore inclined to the plane of the binary by 
34.5°. What determines this preferred inclination will be 
the subject of future investigation of the multiple com- 
ponents of this system. 

Our fit also tells us the angle of rotation of the plane 
of the binary from north-south. Since we assume a sym- 
metry in the precession of the ruff, spending equal time 
rotated clockwise and anticlockwise from the plane of the 
binary, the offset of the fit tells us that the plane of the 
binary (into the sky) is rotated ~ 2° east of north. 

This angle is comparable with, but not equivalent to, 
the mean angle perp endicular to SS433's radio jets of 



Johnston|[l98l| ). 

arped and possibl 



10° east of nort h (Blundell & Bowler 2004 Hjellming 



But since the jets are produced by a 
warped and possibly non-co-planar accretion disk around 
the compact object we think this disagreement may be 
informing us about misalignment of the disk plane and 
jet-launch axis. 



6. 



INDEPENDENT INFORMATION ON PHYSICAL 
PARAMETERS 



6.1. Independent constraints on SS433's orbital 
eccentricity 

Optical photometry, such as in |Kemp et ah] fll986| ), 
reveals that SS433's primary and secondary minima 
eclipses lie close to 0.5 of an orbital phase apart. This 
points towards a circular (non-eccentric) orbit, but is not 
conclusive because the orbit could still be eccentric if the 
axis of symmetry 5 is along our line-of-sight. 

There are other pointers towards a non-zero eccentric- 
ity in SS433's orbit. Fi rst, the sinusoidally varying com- 
ponent of th e jet speed ( |Blundell fc Bow ler 2005 Blundell 
et al.|[2007 ) depends on orbital phase, so there is some- 
thing jEaTbreaks the symmetry of the orbit. Second, 
Perez & Blundelll (120091) find that the radius of the com- 



panion star must be rather large (> 39 R©) if the orbit 
is circular; but smaller radii for the star are permitted if 
the orbit has some eccentricity. 



6.2. Independent confirmation of the mass fraction of 
SS433's binary components 

Blundell et al. (20081), from over a month of high reso- 



lution optical spectroscopy of the circumbinary disk, ex- 
tracted the masses of SS433's binary components as be- 
ing ~ 16 M for the compact object and ~ 24 M for 
the companion. Their measurements were also consistent 
with 18 and 22 M respectively, so on the basis of this 
independent determination a mass fraction of between 
0.4 and 0.45 seems likely. 

We remark that the changing orientation of the pre- 
cessing circumbinary disk will be manife sted in how split 
the c ircumbinary lines appear to be ( |Blundell et al.| 
2008|. When the circumbinary disk is maximally face- 



o n to Earth, then the line-splitting will be reduced (see 
§5.1), and with insufficient spectral resolution, becomes 
harder to resolve. 

7. CONCLUSIONS 

Comparing the results of our simulations (Q with ob- 
servations of SS433's ruff angle (Fig [2j) we conclude that 
to achieve a short ruff precession period of ~ 42 binary 
periods it is necessary that the circumbinary material or- 
bits very close in to the binary, i.e. a radius ~ twice the 
binary semi-major axis. 

In addition, it is necessary to have a high binary mass 
fraction, i.e. binary components of similar mass, and an 
eccentric binary orbit. 

Our simulations show that inclined circumbinary or- 
bits are surprisingly stable throughout binary mass- 
fraction/eccentricity parameter space. The motions of 
collective orbits, comprising circumbinary disks, have 
considerable potential to reveal system properties, if suit- 
ably time-resolved, high-resolution data are available. 
For the case of SS433, its changing ruff angle seems to 
confirm a mass fraction close to equality and a non-zero 
eccentricity. 
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